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I review our proof that long range forces induced by the exchange of massless neutrino- 
antineutrino pairs do not affect the stability of neutron stars. 
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1 Introduction 



Consider two test neutrons at rest, separated by the distance d 3> 1/Mz- If neutrinos are massless, they 
induce (see Fig. [IJ) an effective long range interaction between the two neutrons which, on dimensional 
grounds, must fall like 

(Gp ~ W~ 5 GeV~ 2 is the Fermi coupling.) §. This interaction is very weak: at short distances, 
d < 10~ 8 cm, it is dominated by the direct exchange of a Z boson; at large distances, d > 10 _8 cm, 
the gravitational interaction (!) between the two test neutrons is dominant. Similarly, in the presence 







Figure 1: Exchange of a neutrino-antineutrino pair between two neutrons (n) at rest. 

of N neutrons (taken to be at rest, to simplify), long range, multibody interactions are induced by 
the exchange of a neutrino-antineutrino pair between any subset of k neutrons, k < N. |] Consider 
then the contribution of these interactions to the total potential energy of the N neutrons. As each 
coupling to a neutron brings a factor of Gf, the k-body contribution gets smaller as k increases. 
However, each contributions has to be added with the adequate combinatoric factor, that, at least 
superficially, is proportional to the number of ways to take k neutrons out of N, oc N k for large N. 
To be specific, consider a neutron star, that contains N ~ 10 57 neutrons - corresponding to about 1.4 
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Figure 2: Expansion of the self-energy in k-body neutrino-antineutrino exchange. 

solar mass- within a volume of radius R ~ Wkm, Fig. [3]. The self-energy of the neutron star from 
neutrino exchange is naively represented by the serie of Fig. [2], where the dots now mean insertions 
of the neutron density, pjy. The term with k insertions of the neutron density scales approximately as 

(i) 

where is a dimensionless numerical coefficient. The dimensionless parameter that governs the 
expansion is thus GfN/R 2 ~ GfPuR- For a neutron star, this parameter is O(10 12 ) (!). 

Attempts to a direct summation of this series (truncated after N terms) have yielded enormous 
values for the interaction self-energy §, which led some authors to claim that neutrinos must have a 
mass of at least 0.4 eV in order to allow neutron stars to exist Qtirtlf] Our approach to this problem 

b One-loop diagrams with an odd number of neutron insertions vanish. PI 
c This problem has an interesting story. It was first raised by R. Feynman in his Caltech lectures on gravitycl. He asked 
himself whether long range neutrino forces could mimic gravity, but ruled it out because the equivalence principle was 
not satisfied. He noticed that a large quantity of matteri-iwould render the expansion in k-body interactions divergent. 
Hartle addressed this problem in a cosmological context u. He noticed that the large scale repartition of matter in the 
Universe is such that the perturbative expansion is actually well-defined and leads to totally negligeable effects. He also 
remarked that the long range neutrino interactions have no effect in presence of an homogeneous distribution of matter. 
Fischbach later recognized that the is expansion parameter is S> 1 inside of a neutron staru. 



Figure 3: Exchange of neutrino-antineutrino pairs between the neutrons of a neutron star. 



is different. Following Abada et al, we claim that the series represented in Fig. [2] is meaningless 
if GfPnR <; 1, and must be resummed in a non-perturbative way to get a sensible result i'l. We 
have performed this resummation and find that the apparent infrared divergence is an artifact of the 
expansion and that the energy is finite and well-behaved. This implies in particular that there is 
no lower bound on the mass of the neutrinos. Along the way we have encountered some interesting 
physics. In particular, we show that the ground state of the system has a non-zero neutrino charge - a 
result which was previously anticipated by Loebi. In our simple model for the density of the neutron 
star it is straightforward to calculate both the energy and the neutrino number of a neutron star. Let us 
mention that several groups have examined various aspects of this problemi'00llUll0. In particular, 
Arafune and MimuraEJ have confirmed our asymptotic result using an analytical approximation. 



2 The neutrino ground state in a neutron star 



We want to compute the contribution to the self-energy of the neutron star due to long range neutrino 
forces, or equivalently, the shift in the neutrino ground state energy in presence of a neutron star. 
TheJatter can be defined in terms of the neutrino Hamiltonian iJ( ) i n the presence (absence) of the 
starlit: 

W = (6|F|6) - (0|iJ |0). (2) 

Here |0) denotes the neutrino ground state in presence of the star, while |0) denotes the usual neutrino 
vacuum state. As we have already alluded, the state |0) contains in general a non-zero neutrino 
number (i.e., it is "charged"). Note that the expression in Eq. (||) is a formal, ultraviolet (UV) 
divergent quantity which needs to be renormalized. This renormalization may be done using the usual 
techniques of quantum field theory. 

In order to proceed, it is convenient to introduce a low energy effective Lagramdan for the neutrino 
field. After integrating out all of the other particles in the theory, one obtains EJ' 8 " 
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where ipL = h(l — 75)^! and where 



(3) 



a(x) = G F p r , 



20 eV 



(4) 



is the electroweak potential induced by the finite neutron density (typically p n « 0.4 fm -3 in a neutron 
star). This potential is identical to the one which is usually considered in the well-known Mikheyev- 
Smirnov-Wolfenstein (MSW) effect til The potential term in Eq.(|3|), which is attractive for neutrinos 
and repulsive for antineutrinos, resembles a position-dependent chemical potential, and it is maybe 
not surprising that the ground state of the system can carry a non-zero neutrino number]] ^ 

d This conclusion could be premature. Forr-i« stance, in 1 + 1 dimension, such an external potential has absolutely no 
physical effect (if the neutrinos are massless) Eil. 

e From the effective Lagrangian point of view, it is manifest that an homogeneous distribution of neutrons, p n = const, 
has no net effect: it can simply be gauged away. (See Hartle's comment in footnote 6.) 



The Schwinger-Dyson expansion for W in terms of the potential a(x) and the neutrino propagator 
Go(x, x';uj) gives! 
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This series corresponds precisely to the one which is represented diagrammatically in Fig. [2]. A useful 
aspect of the perturbative expansion is that it isolates the UV divergence in W . In particular, the 
only UV divergent term in Eq. (j6j) is that with k = 2 which is related to the vacuum polarization of 
the Z boson in the complete theory. While the terms with k > 4 are separately UV finite, their sum 
is ill-defined for aR > 1. The non-perturbative resummation of these terms is the main goal of our 
calculation. 

2.1 The 2-body contribution. 

We begin by an estimate of the two-body contribution to the self-energy of a neutron star, W%. From 
Eq.(||) this contribution is UV divergent, W2 ~ A 2 , with A some UV cut-off . The procedure to 
compute this term is in principle well-defined. The diagram can be computed and renormalized both 
in the effective theory and in the underlying theory - i.e. the Standard Model - and the respective 
results must be matched at some scale. A natural scale for matching is the inverse size of the neutrons, 
A ~ lGeV. From the effective Lagrangian point of view, this implies that the coarse grained structure 
of the star is taken into account. At larger distances, d ^> GeV~ l it is however sufficient to consider 
a smooth, continuous distribution of neutrons. This simplification is helpful when we come to the 
evaluation of the higher order, UV finite, multibody contributions. 
A rough estimate for W2 isQ&Q£3 

W 2 ~ + a 2 A 2 R 3 (7) 

The two-body interactions give a contribution oc i? 3 , and proportional to the number of neutron pairs 
in the star, oc p 2 ~ A^ 2 . Also, because the two-body interaction is repulsive, this contribution is 
positive. For a neutron star, this estimate gives a negligeable contribution, 

W 2 ~ 10 16 kg < 1.4 M ~ 10 31 kg. (8) 

2.2 The k-body contributions, k > 4. 

The terms in the Schwinger-Dyson expansion with k > 4 are UV finite. However, there are sensitive 
to the size of the star, oc (a R) k . For aR > 0(1) these terms must be resummed in a non-perturbative 
way. We begin by making some useful simplifications. First, to study long distance, UV finite effects, 
we can forget the neutrons and only consider their mean field effect: a(x) is taken to be a smooth 
function of x. Furthermore, we approximate the shape of the star by a spherical square well potential, 
with depth a ~ 20eU and radius R. For k > 4, the only dimensionless parameter is aR. If aR < 1 the 
Schwinger-Dyson expansion is well-defined (the "weak coupling regime"), while it must be resummed 
if aR > 1 (the "strong coupling regime"). If we want to understand the transition from the "weak 
coupling" to the "strong coupling" regime, we only have consider systems of rather small size: e.g. 
for fixed a ~ 20eV , 0.1 < aR < 100. The extrapolation to R ~ 10km will be trivial. 

To resum the k > 4 terms, we use the following expression -adapted from Schwinger for the 
shift of the neutrino vacuum energy in presence of the external electroweak potential: 

W = — ]T(2Z + 2) / dw [S t (u) + S^-u)] . (9) 
2vr^ Jo 

where 8i(rLu) is the scattering phase shift of an incident neutrino (antineutrino) of energy (—)u and 
I labels the orbital angular momentum. The factor (21 + 2) = (2j + 1) is the degeneracy factor for a 
given energy ±u and total angular momentum j. 



The details of our calculations can be found inl3. Here, I only summarize the main results. 



- Eq. (P) is a formal [TV-divergent expression which need to be renormalized. We achieved this 

by substracting the leading Wi = 0(a 2 ) term in the Born expansion of the phase shifts. The 
"renormalized" term Wi is then like in Eq.(Q). 

- If the neutrino are massless, there are no bound states, only resonances. This is because a neutrino 

must be able to flip chirality i>l — ► vr in order to form a bound state. Resonances exists because 
no chirality flip is necessary to bend the trajectory of a massless neutrino. (See Loebi.) These 
resonances become extremely narrow - i.e. long lived - as R increases. 

- The physical implication of these resonances is that a non-zero neutrino charge can be "confined" 

within the potential/neutron star. Intuitively, the external electroweak potential induced by 
the finite neutron density polarizes the neutrino vacuum. If the external field is strong enough, 
i.e. aR > 1, the neutrino vacuum is unstable toward the creation of neutrino-antineutrino 
pairs. In an open system, like a neutron star, the antineutrinos fly away to spatial infinity while 
the neutrinos effectively "screen" the external electroweak potential. This effect is obviously 
non-perturbative in the external potential. 
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Figure 4: Renormalized exact charge, Q/Q v , and self-energy \(W — W2)/W u \ as a function of aR. 



As R increases, we observe a first order phase transition at aR ~ 1 between the "weak coupling" 
and the "strong coupling" regimes, after which the self-energy rapidly reaches an asymptotic 
behaviour, Fig. [4]. The calculations show that this asymptotic value coincides with the ther- 
modynamic limit of a ideal Fermi gas of massless neutrinos with chemical potential a, trapped 
within a volume V = 4/37ri? 3 : 

a 4 R 3 

W-W 2 « W v = -—- 

187T 

Q « Qu = 2 { ^- (10) 



This behaviour for very large aR was anticipated by Abada et alu and confirmed by Arafune 
and Mimura0. 

- Because there is no other dimensionless parameter in the theory than aR, it is trivial to extrapolate 

these results to a realistic neutron star, R ~ Wkm. There are about 

Qu ~ 10 36 (11) 

low energy, u ~ eV, neutrinos trapped in a neutron star. Expressed in kilograms, they contribute 
about 

W - W 2 ~ -30kg (12) 

(for each generation of massless neutrinos) to the mass of the star, to be compared to a total 
mass of 10 31 A;^. 

- This effect is extremely small, probably of no observable consequence. Incidentally, this proves that 

massless neutrinos can not affect the gravitational stability of a neutron star. 



3 Conclusions 



The physics underlying the resummation problem is non-perturbative. The neutrons are the source 
of a neutral electroweak effective potential that polarizes the neutrino vacuum. When the system is 
large/dense enough (the "strong coupling" regime), like in a neutron star, the energy of the system 
is lowered by spontaneous creation of neutrino-antineutrino pairs. In an open system, again like a 
neutron star, the ground state contains a "neutrino condensate" . These results show that there is no 
"mysterious" long range neutrino force at work in a neutron starEj. 
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